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Abstract 



We consider the direct detection of weakly interacting massive particles 
(WIMPs) reaching the Earth from outside the Milky Way. If these WIMPs 
form a distinct population they will, although of much lower flux than typical 
galactic halo WIMPs, have a number of features which might aid in their 
ultimate detectability: a high and essentially unique velocity (~ 600 km/s in 
the galactic rest frame) due to their acceleration in entering the Milky Way, 
and most likely one or two unique flight directions at the Earth. This high 
velocity may be experimentally advantageous in direct detection experiments, 
since it gives a recoil signal at relatively high energy where background is gen- 
erally much reduced. For a density of extragalactic WIMPs comparable to 
the critical density of the universe the count rate expected is very roughly the 
same as that of fast galactic WIMPs. If there is an increased density rela- 
tive to critical density associated with the Local Group of galaxies, say 10-30 
times the critical density, there is a corresponding increase in rate and the 
extragalactic WIMPs would show up as a high energy shoulder in the recoil 
energy distribution. Evidence of such WIMPs as a separate population with 
these distinct properties would offer interesting information on the formation 
and prehistory of the galaxy. 
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I. INTRODUCTION 



The nature of the dark matter in the universe is one of the outstanding questions in 
astrophysics and particle physics. One of the favored candidates is the WIMP (weakly 
interacting massive particle). These particles are presumed to have masses in the range 1 
GeV/c 2 to a few TeV/c 2 . 

Much work has been done studying the possibilities for detecting these particles. Pos- 
sibilities include direct detection |IJ, where the particle interacts with a nucleus in a low 
temperature detector, and is identified by the keV of energy it deposits in the nucleon; and 
indirect detection, whereby the particles are captured in the Sun or Earth, sink to their 
centers, annihilate with one another in the core, giving rise to particles including neutrinos 
which can be detected by experiments on the surface of the Earth 0,0]. Alternatively, the 
particles annihilate in the galactic halo |4| or the galactic center || and produce anomalous 
components in cosmic rays. 

Most of this work has focused on detection of WIMPs that belong to the virialized halo of 
the Milky Way galaxy. However, some authors [|6| P~0j| have started to consider the detection 
of WIMPs that are not yet virialized with the Milky Way halo, specifically because they have 
fallen into the Milky Way galaxy only relatively recently. The general result is that late- 
falling WIMPs reach the solar neighborhood as streams of cold dark matter particles with 
small velocity dispersion and high bulk velocity. Only few of these streams are expected, 
each one coming from a particular direction. 

In this paper we discuss the possibility of detecting "extragalactic WIMPs." These 
would appear as streams of WIMPs passing the solar neighborhood. The existence of a 
separate dark matter population with properties distinct from that of galactic WIMPs is 
plausible although not certain. Its presence depends on the prehistory and the evolution of 
the galaxy. So verification of such an extragalactic component would provide very interesting 
information on the history and dynamics of the galaxy. 

We address the question of the detectability of extragalactic WIMPs within a simple 
model in which the galaxy is moving through a background field of ambient WIMPs. In this 
way, there is a simple distinction between bound and extragalactic WIMPs. Although this 
is not the commonly used picture at the present time - where the galaxy is still in a state of 
formation and there is no clear distinction between bound and extragalactic WIMPs - our 
infalling ambient WIMPs may resemble the streams found in detailed simulations. Indeed, 
the investigation of our simplified model provides the opportunities to look at possibilities 
for detection and to begin to investigate this complex of questions. 

Any such extragalactic WIMPs reaching the Earth from outside the halo will have been 
accelerated considerably by falling in the potential well of the Milky Way. They may thus 
be expected to have the local escape velocity from the Galaxy, roughly 600 km/s, to be 
compared with the 200-300 km/s typical of halo WIMPs. Furthermore, in a picture where 
the galaxy is moving through a background of WIMPs essentially at rest, they will appear 
in the rest frame of the galaxy as a "wind" with definite incoming direction and velocity 
at infinity. In the case of a spherically symmetric halo, this leads to their appearance at 
the Earth from one or two essentially unique directions, one direction corresponding to 
WIMPs on their way into the galaxy, the other to WIMPs on their way out of the galaxy. 
Hence in this picture one has extragalactic WIMPs at the Earth with the distinct features 
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of essentially unique energy and direction of motion. 

In direct detection based on the observation of nuclear recoil, the unique energy has two 
interesting implications. First, the recoil energy for extragalactic WIMPs is higher than that 
for typical halo WIMPs. Since, in general, the background decreases sharply with energy 
|TT| , this may be a very favorable point experimentally. Secondly, a monoenergetic incoming 
flux produces a characteristic recoil spectrum quite different from the recoil spectrum from 
galactic WIMPs. For example, if the WIMP is not too heavy with respect to the target 
nucleus, the recoil energy spectrum has a "box-like" shape, quite different from the tail of a 
Maxwellian distribution. 

Concerning the definite arrival direction at the Earth, exploitation of this should even- 
tually be possible with direction-sensitive detectors Hl2| . The directional sensitivity is of 



particular interest in connection with extragalactic WIMPs, and indeed could be used as a 
filter to increase the extragalactic WIMP signal relative to halo WIMPs and background. 

When directional sensitivity becomes available flE^ , an improvement of the signal can 
be hoped for, in addition to interesting additional information on the extragalactic WIMPs 
population. Unfortunately with present knowledge it is not possible to calculate this di- 
rection at the Earth, since it depends on the unknown shape of the Galactic gravitational 
potential (see below). On the other hand were a definite direction ever to be verified, this 
direction, like the feature of a distinct energy, would provide very interesting information 
on the structure and evolution of the galaxy. 

Naturally the absolute rate to be expected from extragalactic WIMPs is very small and 
this is undoubtedly the major obstacle to a practical realization of these ideas. However with 
progress in background reduction and increasing mass of cryogenic detectors it is perhaps of 
interest to consider detection of extragalactic WIMPs for the future in a preliminary way. 
We shall find that under the assumption of an extragalactic WIMP population density of 
the order of the universal critical density, the expected count rate is very roughly the same 
as that of fast halo WIMPs. If there is an enhanced density associated with the Local 
Group of galaxies, say by a factor of 10-30 (see below), then the rate is correspondingly 
increased. Thus a distinct monoenergetic WIMP population -if it exists- could show up as 
a shoulder above the halo WIMP spectrum. Verification of such a component would provide 
very interesting information on the history and dynamics of the galaxy. 

We proceed as follows. First, in Section O, we discuss the nature of the incoming flux 
of extragalactic WIMPs for an observer in the Solar System, i.e., the velocity, the direction, 
and the flux enhancement due to gravitational focusing. Then, in Section [TTT|, we discuss 
interactions in a detector: cross sections and recoil energy spectrum. In Section [IV], we 
discuss the expected density of fast galactic WIMPs, and how we will use it as a basis for 
comparison. In Section |V|, we discuss the expected count rates due to our galaxy moving 
through a critical density of extragalactic WIMPs. Finally, in Section VI, we discuss the 
expected count rate due to Local Group WIMPs overdense by a factor of 10-30 relative to 
critical. We present our numbers relative to the expected count rates for fast halo WIMPs; 
the precise count rates depend on the details of the detectors, details which are the same 
for Halo and extragalactic WIMPs. 

A general note for this paper: given the level of uncertainty tolerated, we neglect the 
small velocity of the Earth in its motion around the Sun, and also neglect the peculiar 
velocity of the Sun with respect to a frame rotating with the galaxy at the Sun location (the 
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so-called local standard of rest). In other words, we take the Earth and the Sun as moving 
at the rotation velocity of the galaxy at the Sun's position, v Tot = 220 km/s. With these 
considerations in mind, we will refer to Earth, Sun, and Solar System interchangeably. This 
simplifying assumption precludes us from considering the annual modulation of the signal 
due to the Earth motion around the Sun [[n||14]]. We only comment that the amplitude of 
this modulation is approximately equal to the ratio of the Earth's speed, 30 km/s, to the 
WIMP speed, and so decreases as the speed of the WIMPs increases. Hence for extragalactic 
WIMPs arriving on Earth at the escape velocity, the amplitude of the modulation is expected 
to be approximately 5%. 



II. LOCAL FLUX OF EXTRAGALACTIC WIMPS 

In this section, we discuss the flux of extragalactic WIMPs as seen by an observer moving 
with the Solar System. The WIMP flux is an essential quantity for estimating the total 
interaction rate, which is proportional to the flux F and the interaction cross section a, 

Rtot = Fa, (1) 

where F = nw is the product of WIMP number density times velocity. However, the 
WIMP speed itself is also significant in that it governs the recoil energy in the detector, 



and so determines the differential interaction rate (see Section [HI] ). Furthermore, given the 
prospects for a directional WIMP detector [15], also the direction of the WIMP velocity is 
of interest. 

Velocity: When an extragalactic WIMP enters the Milky Way, it is accelerated to at 
least the local escape velocity f esc by the potential well of our galaxy. The local escape 
velocity from the Milky Way at the position of the Sun is not well determined. Following 
the discussion in |TJj], we take a value of about t> esc = 600km/s. The local speed v of an 
extragalactic WIMP with velocity Voo far away from the Milky Way follows from energy 
conservation 



+ v 2 csc . (2) 

Since we anticipate that WIMPs outside the galaxy have rather low velocities compared to 
the escape speed, to a first approximation WIMPs of extragalactic origin at the Earth will 
have a unique velocity of v ~ t> esc in the rest frame of the galaxy. In a simple model we 
consider below, the increase in v over t> esc due to an initial velocity of the WIMPs is on the 
order of a few percent. 

The transformation from the WIMP velocity in the galactic rest frame to the WIMP 
velocity relative to the Earth is obtained by a simple vectorial addition of its velocity relative 
to the galaxy and the galactic rotation velocity. As such, it depends on the arrival direction 
of the WIMP when it approaches the Solar System. We define tp to be the angle between 
the following two vectors: the vector starting at the Earth and pointing in the direction the 
WIMPs are moving near the Earth, and the vector starting at the Earth and pointing in 
the direction of the galactic rotation velocity. Then the WIMP speed relative to the Sun is 



w = Jv 2 + v? ot - 2vv rot cos if). (3) 
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Arrival Direction: Statements on the arrival direction of extragalactic WIMPs are model 
dependent. We will focus on a plausible model in which the Milky Way is in motion with 
respect to the extragalactic WIMPs in its vicinity. Viewed in the rest frame of the Milky 
Way, there is then a parallel flux of extragalactic WIMPs entering the galactic halo with a 
definite velocity at infinity Voq. (Section [VT] discusses a particular model of this type in more 
detail.) 

The further progress of the particles, and in particular their direction upon reaching 
the Sun, depends very much on the details of the potential in the galactic halo. In a 
flattened halo, WIMPs may follow complicated non-planar trajectories. We will consider 
instead a spherical potential. In this case, conservation of energy and angular momentum 
impose that the WIMP trajectories lie on the plane defined by the Sun, the direction of 
galactic motion, and the galactic center. Spherical symmetry would also impose that only 
a discrete number of trajectories on the orbital plane intersect the location of the Earth. 
Thus the detector should see WIMPs coming from a few unique directions, which would 
be somewhat blurred by the velocity dispersion of the extragalactic WIMPs. On the one 
hand, it is not easy to predict what these directions are, because they strongly depend on 
the radial behavior of the galactic gravitational potential. (Some examples will be given 
later in Section [VT| and in the Appendix.) On the other hand, detecting and determining 
the direction of possible extragalactic WIMPs would provide very interesting information 
on the gravitational potential of the galactic halo. 

Focusing: A further interesting effect of the galactic potential is gravitational focusing, 
which will tend to increase the flux at the Sun relative to its value at infinity. Again, specific 
results even in the simple model depend on the details of the galactic potential. 

However it is interesting that one may find a simple result for a kind of average enhance- 
ment by the following argument. We view the galaxy, in its rest frame, as being subject to 
a parallel flux of particles, which at great distances is Foo = PooVoo/m, where p^ = p ex tragai 
is the density of extragalactic WIMPs far away from the Milky Way and m is the WIMP 
mass. The flux is the appropriate quantity to study since it is divergenceless in a stationary 
situation. Consider a particle incident at very great impact parameter. In its trajectory 
through the galaxy, or rather the galactic halo, it will remain far from the center and suffer 
a very small deflection. As we reduce the impact parameter, the trajectory will pass closer 
and closer to the center. Consider the sphere around the center of the galaxy on which the 
Sun is located, of radius -Rsun- As the incident impact parameter is reduced, we will reach a 
critical value of the impact parameter where the trajectory of the particle just touches the 
sphere of radius -Rs un - 

Now every incident particle with less than this impact parameter, which we call & max , will 
pass through the sphere. Hence the total number of particles per second passing through 
the entire surface of the sphere is 



If we divide this by the surface area of the sphere we will obtain the average flux over 
the sphere: 



sphere 




(4) 
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We now use angular momentum conservation to find b mSLX . The angular momentum per 
unit mass as calculated at infinity for the particle on the critical trajectory is Woo^max- On 
the other hand when it just touches the sphere it is vRs nn . 

We thus obtain 



and the average flux 



v 

&max = -Rsun (6) 

Voo 



V 2 

-^sphere ^'oc 1T~ • (7) 



This equation is invalid when Voo = 0, in which case there is no WIMP wind and the simpli- 
fied model described here should be amended by including the WIMP velocity dispersion.) 

We can also express this result in terms of WIMP densities instead of fluxes. From 
-^sphere = PsphereV / m and = p^v^/m, it follows that the WIMP density averaged over 
the sphere of radius i?sun is given by 

v . . 

Psphere Poo [y ) 

Voo 



With v ~ f esc ~ 600 km/s and ~ 75 km/s (see Section |VI|), we obtain a magnification 

factor p S ph e re/Poo ~ 8. 

This is an essentially exact result within the model but unfortunately it is only the av- 
erage over the entire sphere in question, while instead we would like to know the density 
of extragalactic WIMPs at the location of the Solar System p ex tragai,iocai- This density de- 
pends on the details of the gravitational potential, and is in general correlated with the 
arrival direction of the extragalactic WIMPs at the Sun. Obtaining the exact value of the 
enhancement factor at the position of the Earth is quite complicated. Given the greater 
uncertainties of other numbers here, this calculation is not currently justified. We therefore 
parametrize the uncertainty in the focusing enhancement by a factor 6f oc , and write 

Pextragal. local ^focPsphere- (9) 

Since the position of the Sun on this sphere is not exceptional^ the average should roughly 
characterize the enhancement factor to be expected. In more detailed calculations using 
various potential shapes, we find there can be variations of two or three in the enhancement 
factor but, barring very special conditions, not of an order of magnitude; i.e., we find 
0.3 < b ioc < 3. 



1 As discussed in Section VI, our best guess is that WIMPs enter the Galaxy from the direction 



of Andromeda, which lies ~ 120° from the Galactic center as seen from the Solar System. 
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III. INTERACTION RATE IN A DETECTOR 



The total interaction rate of WIMPs with a target detector nucleus is given by 

Rtot = Fa, (10) 

where F = pw/m is the WIMP flux, m is the WIMP mass, and a is the cross section for the 
interaction. Since the flux of Milky Way WIMPs largely dominates over any imaginable flux 
of extragalactic WIMPs, the total interaction rate in a detector is expected to be largely 
dominated by galactic, instead of extragalactic, WIMPs. 

However, as seen in Section f|, the local speed of an extragalactic WIMP, and so its 
kinetic energy, is higher than that of any galactic WIMP. In scattering with a detector 
nucleus, the more energetic extragalactic WIMP can deposit more energy, and so impart a 
higher recoil energy to the detector nucleus than any galactic WIMP. Hence selecting events 
with high recoil energy helps in searching for a signal due to extragalactic WIMPs. It is 
therefore important to consider the interaction rate differential in the nucleus recoil energy, 
and compare this differential rate for extragalactic WIMPs with that expected from fast 
galactic WIMPs. In the rest of the Section, we review the generalities of the differential 
interaction rate, that is the recoil spectrum, and in the following sections we discuss and 
compare the recoil spectra expected from extragalactic and fast galactic WIMPs. 

In general, the differential cross section for elastic scattering of a WIMP of mass m with 
a nucleus of mass M can be written 



where w is the speed of the WIMP relative to the detector, q is the momentum transfer 
in the elastic collision, fi = mM / (m + M) is the WIMP-nucleus reduced mass, and J-"(q) 
is a nuclear form factor that takes into account the loss of coherence in WIMP-nucleus 
interactions for momentum transfers comparable to or larger than the inverse nuclear radius 
(we normalized .F(O) = 1). 



For purely scalar interactions, 



4/i 



2 



^.scalar = — [Z f p + (A - Z) f n ] 2 . (12) 



7T 



Here Z is the number of protons, A — Z is the number of neutrons, and f p and f n are the 
neutralino couplings to nucleons. For purely spin-dependent interactions, 

<r ,spin = (32/ 7r)G 2 F fi 2 A 2 J (J + 1) . (13) 

Here J is the total angular momentum of the nucleus and A is determined by the expectation 
value of the spin content of the nucleus (see [|ll p!3i^5H T7|]). 



The nucleus recoil energy equals the energy lost by the WIMP and is simply 

E=-^-. (14) 
2M V ; 
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The maximum momentum transfer is g max = 2/iv, and the maximum recoil energy imparted 
by a WIMP of speed w follows as 



as 



The recoil spectrum for WIMPs of a given speed w then follows from a change of variables 
dR da \ P °° F{q), for E < ff max , 

— = 2F— = I 2fi 2 mw (16) 
dE ^ U for E > E mSLX . 

We have used F = pw/m, and we have divided by the nucleus mass so as to obtain a rate 
per unit detector mass, as customary 

As argued in Section |TJ, extragalactic WIMPs in the solar neighborhood have essentially 
a single speed w. Inserting the focusing enhancement factor as in eq. (§), we obtain that 
extragalactic WIMPs give rise to the following recoil spectrum 



dR 

Je 



hoc V PooCXo , , „ 

7T^— F{Q), for E < ^max, /, 7 n 

w Voo 2/j, z m yi-i) 

extragal I 0, for £ > E max . 



Here w is given by the expression in eq. (Q). Notice that both the focusing enhancement 
factor bf oc and the WIMP speed relative to the Sun w depend on the details of the galactic 
gravitational potential. 

The nuclear form factor T{q) depends on the type of WIMP-nucleus interaction and 
on the mass and spin distributions within the nucleus. Some of these distributions can be 
found in |jl5| -|P7jl. For the present simple considerations, we consider only the case where 
form factor effects are negligible, q <ti h/R, where R is the radius of the nucleus. (For a 
velocity of 600 km/s this translates into either m < M/[(Mc 2 /20GeV) 4/3 - 1] or M < 20 
GeV/c 2 independently of m.) The cross section is then simply a . 

With extragalactic WIMPs all moving at a single velocity and form factor effects ne- 
glected, the differential cross section is constant up to the kinematic maximum. Thus from 
eqs. ( PD and flUD, for low mass extragalactic WIMPs, the recoil spectrum in a detector is 
"box-like", i.e. dR/dE = (pao)/(2p 2 mw) = constant is flat as a function of recoil energy E 
up to the maximum recoil energy E max given in eq. (|I5|) . This signature is completely differ- 
ent from that obtained from a Maxwellian velocity distribution; for halo WIMPs, one must 
average these boxes over different velocities. We recall the procedure in the next Section. 



IV. COMPARISON WITH FAST GALACTIC WIMPS 

It is useful to have a simple comparison of the expected detection rates for extragalactic 
and galactic WIMPs. For galactic WIMPs, we use the detection rates previously calculated 
by many authors as a basis for the comparison. In more detail, we compare the extragalactic 
rates with the galactic rates obtained under the simple assumption that the galactic WIMP 
speeds obey a Maxwellian distribution in the rest frame of the galaxy, without truncating it 
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at the escape velocity. This is the WIMP velocity distribution used in the earliest analyses 
of direct detection by Drukier, Freese, and Spergel |13| and Freese, Frieman, and Gould [14 
We write it as 



fMw( v ) dv = 47r PMW 



1 3/2 



2lTV 



1 1- 



v exp 



3v 2 



2v 



1 1- 



dv , 



(18) 



where we take the dispersion velocity v = 270 km/s. Assuming an isothermal sphere, the 
local galactic WIMPs have a density of ~ 0.3 GeV/cm 3 /c 2 , which amounts to 



Pmw ~ 60,000 p crit 



(19) 



(for a Hubble constant of 70 km s _1 Mpc -1 ). Here p cr it = 5 x 10 _6 /?-7 GeV/cm 3 /c 2 where 
hjQ is the Hubble constant in units of 70 km/s/Mpc. 

The distribution in eq. ( p~8| ) that we use as a basis for comparison, although useful for 
this purpose, may not be the real distribution of WIMP velocities, first, because it allows 
galactic WIMPs with speeds in excess of the escape speed, and second, because the real 
nature of the WIMP velocity distribution in the galaxy is unknown. The effects of other 
halo distributions on direct detection have been studied for example in refs. |7 HT0| . [T8 . 



We now present two ways to estimate the ratio of extragalactic to galactic event rates. 
One method gives an order of magnitude estimate, the other is more detailed, but they give 
essentially the same result. 

The quick method is based on the comparison of the galactic and extragalactic WIMP 
fluxes at the escape speed in the galactic rest frame, f esc ~ 600 km/s. How many galactic 



WIMPs are moving this fast? We use the Maxwellian velocity distribution of eq. ([18]) to 
find that the fraction of fast galactic WIMPs (out of the total galactic WIMPs) that are 
within Av = 1 km/s of escape velocity (corresponding to 1% energy resolution) is given by 



Pmw,. 



fi 



so that 



MWl 



Pmw,< 



)Av ~ 4 x 10" 



3 Pcrit- 



Pmw 



(20) 



(21) 



With a reasonable guess of 3 for the focusing enhancement factor, the local density of 
extragalactic WIMPs is 



Pextragal, local ~ 3 p c 



Hence 



Pextragal, local ~ PMW,esc 



Pcrit 



(22) 



(23) 



The ratio of interaction rates for extragalactic and galactic WIMPs moving with escape 
speed then follows as 



-^extragal Pc 



-RMW,esc Pcrit 



(24) 
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We remind the reader that, in the above equation, we have used the integrated rates for 
extragalactic and galactic WIMPs, with the proviso that for the galactic WIMPs the rate 
is only integrated over velocities near local escape speed (and hence is not the total rate 
integrated over all energies). 

We now discuss the more detailed method of comparing the galactic and extragalactic 
rates. Here we calculate the ratio of differential rates rather than of the integrated rates as 
given in eq. 123. First, we convert the speed distribution in the galactic rest frame to a speed 



distribution as seen by an observer moving with the Solar System. We do this by means of a 
Galilean transformation of the velocities. Denoting with w the WIMP speed relative to the 
Sun, and remembering that we identify the velocity of the Sun with the velocity of galactic 
rotation t> rot , we find 



2nv 



71 2 



1/2 



w 



sinh 



^rot 



3wv 



rot 



V 2 



exp 



3(w 2 + v 



rot J 



2v 



> 1 — 



dw. 



(25) 



The differential interaction rate is then obtained by integrating eq. (16|) over the WIMP 
speed distribution f(w). In general, 



dR = po-p . , 
dE ~ 2fx 2 m W 



fH 



dw, 



W E 



w 



(26) 



where w E = (EM/2/j?) 1 ^ 2 is the minimum velocity of the WIMP required to deposit an 
amount of energy E. For the Maxwellian we use as comparison, we find 



dR 
dE 



MW 



Pmw°~o 
2fi 2 m 



Hi) 



2^rot 



erf 



V3(W E + Vrot) 

V2v 



— erf 



V3(w E - v TOt ) 
V2v 



(27) 



The ratio between the extragalactic and galactic recoil spectra dR/dE\ extraga \ and 
dR/dE\ MW depends on the recoil energy and on the arrival direction of the extragalac- 
tic WIMPs. For our comparison, we adopt the maximum value of this ratio of differential 

rates, 



max 



dR/ (iElextragal 

dR/ dE\wiw 



(28) 



Clearly the ratio of differential rates is maximal when the recoil energy equals the max- 
imum recoil energy E ma _ x . Furthermore, if we neglect the dependence of the focusing factor 
on the arrival direction, the local extragalactic flux is maximal when the WIMP velocity is 
opposite to the direction of galactic rotation, that is ^ = 7r in the notation of Section ||. 
This is therefore the value of TZ max that we adopt for our comparison. Explicitly, we use 



K 1 



where w„ 



v 



Poo 2f ro t 



-'foe 



-'rot 



Vco PMW w r _ 



+ v and v 



erf 



\/3Kax + Vrot) 

V2v 



erf 



V3( 



tUm&x ^rot I 



V2i 



(29) 



uJt, + v 2 Kr . For Voo ~ 75 km/s (see Section |VJ), and with the 

ratio 



values adopted for escape speed, rotation velocity, and velocity dispersion, the 1Z 1 
becomes 
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and using eq. (PL 



TZ max ~ 40,000 b ioc (30) 



1Z max ~0.76 foc -^. (31) 

Pcrit 



This is the ratio that will appear in Figure 1 below. This ratio of differential rates is within 
a small factor equal to the analogous expression for the ratio of integrated rates in eq. (|24"D 
(integrated only over velocities near the escape speed). 

We conclude that within an order of magnitude the ratio of extragalactic to galactic 
WIMP-detection rates is 

rate (extragalactic) 

(32) 



rate (galactic near t> esc ) p crit 



A few comments are in order. The detection rate of galactic WIMPs used for compari- 
son in eqs. (|2~3]), ( |3"T] ) and ( |3"2"D is a fictitious extrapolation of the detection rate of galactic 
WIMPs to WIMP velocities exceeding the escape speed. By definition, no galactic WIMPs, 
i.e. no WIMPs bound to the galaxy, exist with speeds above the escape speed. The com- 
parison ratios in eqs. ([24]) , ( |3~1~D and ([32]) should then be understood not as actual ratios of 
detection rates at specific recoil energies, but only as an informative way to gauge the rates 
due to extragalactic WIMPs with respect to the galactic signal. Indeed, the ratio of extra- 



galactic to galactic WIMP-detection rates may in reality be much larger than what eq. ( p2|) 
indicates: The galactic WIMPs are limited to speeds smaller than the escape velocity, while 
the extragalactic WIMPs may have speeds higher than the escape velocity. Hence the upper 
edge of the "box-like" spectrum of extragalactic WIMPs may extend beyond the highest 
recoil energy from galactic WIMPs. This may make the extragalactic signal stick out above 
the galactic signal more than the numbers in this section (and below) would indicate. On 
the other hand, the ratio of extragalactic to galactic rates may be much smaller than what 
eq. (|32|) indicates: The shape of the galactic potential or the motion of the galaxy with 
respect to the extragalactic WIMPs may conspire to make the extragalactic WIMPs move 
in a direction almost parallel to the direction of galactic rotation at the location of the Sun 
in the galaxy. In this case, the speed of extragalactic WIMPs relative to the Sun would 
be a low 380 km/s, and the recoil spectrum of extragalactic WIMPs would be much lower 
than the recoil spectrum of galactic WIMPs even at the highest recoil energy. This is an 
unfortunate possibility than should be kept in mind. Hopefully, a detector with directional 
sensitivity may be able to help in this situation. 



V. EXTRAGALACTIC WIMPS WITH CRITICAL DENSITY 

We consider smoothly distributed extragalactic WIMPs with density equal to the critical 
density p crit . We expect that the typical velocities of these WIMPs away from our Galaxy are 
small compared to the local escape velocity from the Milky Way. Applying the considerations 
in the previous sections, we find that the density of extragalactic WIMPs near the Earth is 
comparable to the density of galactic WIMPs moving at close to escape velocity (600 km/s), 
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Perifocal ~ PMW.600- 



(33) 



The subscript (crit, local) is intended to remind the reader that this is the local density of 
extragalactic WIMPs that at infinity have a critical density. We find a similar relation for 
the detection rates, 

-Rcrit ~ -Rm\V,600; (34) 

that is, the expected extragalactic WIMP count rate is comparable to that for fast halo 
WIMPs moving at close to escape velocity. 

However, these estimates are very uncertain, because, not knowing the relative motion 
of the galaxy with respect to critical-density extragalactic WIMPs, the considerations at the 
end of the Section |TV| unfortunately apply fully. A somewhat less uncertain case is discussed 
next. 

VI. LOCAL GROUP WIMPS 

It would be interesting if there would be an overdensity (above critical) of extragalactic 
WIMPs near the Milky Way. If so, these extragalactic WIMPs would have a count rate in 
excess of the count rate from halo WIMPs at escape velocity. For example, our galaxy may 
be sitting in an overdensity of WIMPs due to the fact that we are in the Local Group of 
galaxies. If the Local Group would be overdense relative to critical by a factor of 10-30, then 
the expected count rate in a detector would be 10-30 times that of fast halo WIMPs near 
the escape velocity. In addition, if the Milky Way is moving in the direction of Andromeda 
(as argued below), then we would be moving into a wind of Local Group WIMPs. In this 
case, the Local Group WIMPs would approach our detectors from one or at most two unique 
directions. Here we consider this possibility. 

The Local Group of galaxies is a group of galaxies dominated by the Milky Way and 
M31 (a.k.a. the Andromeda galaxy). Currently the number of galaxies stands at 35. For 
a discussion of Local Group properties, see |T9| , pO"| . Each of these big galaxies has a few 
small satellites (these include the Large and Small Magellanic Clouds for the Milky Way). 
In addition, there are some isolated small galaxies. Most of the mass of the Local Group is 
concentrated in the Milky Way and Andromeda. The Milky Way is roughly 400 kpc from 
the center of gravity of the Local Group. The zero-velocity surface, which separates the 
Local Group from the field that is expanding with the Hubble flow, has a radius 1.2-1.8 Mpc 
from the center of gravity of the Local Group. 

The question is, are there any Local Group WIMPs that are not virialized with the Milky 
Way? We will stipulate the existence of a 1O 12 M Local Group WIMP halo (not virialized), 
with the caveat that in the worst case there may be no such halo at all. Recent studies 
of galaxies [f2lH 23"ll indicate that galaxies are larger than previously believed, with radii at 
least 250-300 kpc, and mass-to-light ratios of 145-200. Hence it is possible that much of 
the dark matter in the universe is concentrated in the galaxies themselves. However, even 
if this is the case, it is plausible to believe that when galaxies come into contact with one 
another some of the dark matter is stripped off. Hence it is possible that there are Local 
Group WIMPs not associated with either the Milky Way or M31; these form a Local Group 
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halo but are not virialized. In particular, explaining the rapid motions of galaxies near the 
outskirts of the Local Group may require the existence of a Local Group halo not associated 
with either the Milky Way or Andromeda galaxy. Hence in the remainder of this section we 
will stipulate the existence of a 10 12 M Q Local Group halo of WIMPs. 

We now specialize the considerations of the previous Sections to Local Group WIMPs. 

WIMP velocity and arrival direction: The velocity dispersion within the Local Group is 
o T = 61 ± 8 km s" 1 [[HJ. The Milky Way is moving towards the center of gravity of the 
Local Group with a similar velocity (about 75 km/s). The tangential velocity of the Milky 
Way towards the Local Group has not been measured. However, it is likely to be small: 
the Local Group is an overdense region that was expanding and now collapsing, so that the 
Milky Way is likely to be falling into M31. Hence we can assume the Milky Way is moving 
in the direction of M31 with a speed of 75 km/s. This gives rise to a WIMP wind from the 
direction of M31 with a velocity at infinity ~ 75 km/s. 

The Local Group dispersion velocity of 60 km/s is only 10 percent of the Milky Way 
escape velocity of 600 km/s, and the associated energy is only 1%. Hence energy conservation 
tells us we need to modify the velocity near Earth by only 1%. Thus as above, we find that, 
with respect to the galactic rest frame, the speed of Local Group WIMPs in the solar 
neighborhood is to a good approximation ~ 600 km/s, that is the escape speed from the 
galaxy. 

As discussed in Section ||, the specific arrival direction of extragalactic WIMPs at the 
location of the Sun and Earth depends on the details of the galactic gravitational potential. 
For the sake of illustration, we have assumed a spherical Milky Way halo and we have 
calculated WIMP trajectories for two radial dependences of the gravitational potential: a 
Navarro- Frenk- White potential and a logarithmic potential. (See the appendix for details.) 
We have found that at most two arrival directions at the Sun are possible, one for WIMPs 
on their way into the Milky Way and one for WIMPs on their way out. Thus there are one 
or at most two values of the WIMP speed at Earth. 

Without any constraint on the direction and speed of the relative motion of the Milky 
Way and the Local Group WIMPs, these directions could cover the whole possible range 
from opposite to parallel to the Sun motion in the galaxy (which is roughly the direction of 
galactic rotation). In this case, the WIMP speed relative to the Earth, obtained subtracting 
the Sun velocity from the WIMP velocity vectorially, could have a value (or at most two 
values) in the range 380-820 km/s, which is 600±220 km/s, according to the arrival direction 
at the Sun. 

However, if we assume that the WIMPs are approaching the galaxy at infinity from the 
direction of M31, we find that the WIMP speed relative to the Earth must be larger than 
~ 640 km/s. This gives a higher ratio of Local Group and galactic WIMPs counting rates. 
The reason for the selection of the higher range in speeds is the conservation of the sign of 
angular momentum. To keep track of the angular momentum, it is helpful to introduce an 
observer at the Earth that faces the galactic center with the north galactic pole above his/her 
head. Then M31 lies in the southern hemisphere on his/her left-hand side, at ~ 120° from 
the galactic center (M31 has galactic longitude I = 121?2 and galactic latitude b = — 21?6). 
Now consider the WIMPs that come from M31 and that reach the Earth on their way into 
the galaxy, that is before they reach the point of closest approach to the galactic center. Their 
initial angular momentum points toward the north galactic pole. When these WIMPs arrive 
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at the solar neighborhood, their angular momentum still points toward the north galactic 
pole (conservation of the sign of angular momentum). So for our hypothetical observer they 
possess a velocity component (tangential velocity) directed toward the right, together with 
another velocity component (radial velocity) directed toward the galactic center. The key 
point is that the direction of galactic rotation, and also the direction of the solar motion, 
is to the left of the observer. Hence the angle between the WIMP velocity (sum of radial 
and tangential components) at the solar neighborhood and the direction of galactic rotation 
cannot be smaller than 90°. It follows that the WIMP speed at the Earth cannot be smaller 
than (600 2 + 220 2 ) 1 / 2 = 640 km/s. A more detailed analysis of the geometry also shows that 
the angle between the WIMP velocity at the Earth and the direction of galactic rotation 
cannot be larger than 155°, which gives a maximum WIMP speed relative to the Earth of 
805 km/s. 

Density at the Earth: Using the numbers given above for the total mass 1O 12 M in 
WIMPs and a Local Group radius of 1.2-1.8 Mpc, we find an average Local Group WIMP 
mass density of 

Plgw = (1 - 3) x 10- 28 gm/cm 3 . (35) 

For a Hubble constant of 70 km s -1 Mpc -1 , this amounts to (10-30)p crit . Our galaxy is 
roughly 400 kpc from the center of gravity of the Local Group, about 1/3 of the way out 
to the zero- velocity surface; for a 1/r 2 density profile, we should thus be at roughly 10p cr i t . 
(Since the system is not virialized we cannot be sure of the density profile, but this estimate 
still gives us an indication that is not unreasonable to use an average WIMP density in the 
halo as its value at the outskirts of the Milky Way.) 

Thus the density ratio of nearby Local Group WIMPs to galactic WIMPs moving with 
600 km/s is given by 

PLGW,local ~ (10 — 30)pMW,600 , (36) 

where pLGvyiocai is the density of Local Group WIMPs near the Earth and Pmw,600 is the 
density of Milky Way WIMPs moving with 600 km/s. 

The count rate of Local Group WIMPs in detectors would also be 10-30 times that of 
galactic WIMPs 

-Rlgw ~ (10 — 30)i?MW,600- (37) 

These Local Group WIMPs would have a unique signature in that (i) they come only from 
one or at most two specific directions and (ii) the lighter ones would have a flat recoil 
spectrum. 

Unfortunately the numbers of events in current detectors would still be tiny. The ex- 
pected count rate is 0.1 below that for typical halo WIMPs moving ~ 220 km/s, and 10~ 3 
below the signal from the sum total of halo WIMPs. 

However, as discussed in the Introduction and at the beginning of Section [T| selecting 
events with high recoil energy should help in the detection of Local Group WIMPs, because 
the recoil spectrum of extragalactic WIMPs would show up as a high-energy shoulder which 
may stick out above the spectrum of galactic WIMPs at high recoil energies. We calculate 
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the recoil spectra for extragalactic WIMPs using eq. (|17D and compare it with the galactic 
spectrum from eq. (|27|). As discussed in Section |ITT|, the spectrum of extragalactic WIMPs 
is boxlike. 

Figures 1 and 2 show examples of the expected recoil spectra for Milky Way and Local 
Group WIMPs. For the sake of illustration, we have fixed the mass of the target nucleus 
M = 73 GeV/c 2 , the WIMP mass m = 60 GeV/c 2 , the WIMP-nucleus scattering cross 
section o~o = 10 -35 cm 2 , and we have set the nuclear form factor T(q) = 1. For the galactic 
WIMPs, we have obtained the recoil spectrum from eq. (^) for a Maxwellian velocity 
distribution as in eq. flTg| ) but truncated at the escape velocity (we do not display the 
lengthy formulas). For comparison, the dotted line shows the spectrum due to Milky Way 
WIMPs for the non-truncated velocity distribution we use as a comparison in Section [TV], 
that is eq. (j2T|). For the Local Group WIMPs we have obtained the rate using eq. (|T7|), taking 
frfoc = 1; Plgw — 20p crit . Since we do not know the arrival direction of the Local Group 
WIMPs at the Solar System, we show two possible spectra for these WIMPs, Figure 1 and 
Figure 2. We believe that these two spectra bracket the plausible range. They correspond 
to the two values 90° and 155° for the angle ip between the velocity of the Solar System in 
the galaxy and the velocity of the Local Group WIMPs when they reach the Solar System 
(both velocity vectors starting at the earth). 

In the most favorable case of ip = 155°, the high energy recoils due to Local Group 
WIMPs can be clearly distinguished from the recoils due to galactic WIMPs. In the less 
favorable case of ip = 90° the distinction is more difficult. 

VII. CONCLUSIONS 

Extragalactic WIMPs have several unique signatures which may aid in their ultimate 
detectability. As these WIMPs enter the Milky Way, they are sped by the gravitational 
potential of the Galaxy. Near the solar neighborhood, they are moving with the local escape 
velocity of 600 km/s relative to the galactic rest frame. The relatively high velocity of the 
extragalactic WIMPs can have the advantage of moving the recoil signal in direct detection 
to higher energies, where background is much lower. In addition, the extragalactic WIMPs 
come from only one or two directions. Although it is currently not easy to predict what 
these directions are, because they strongly depend on the radial behavior of the Galactic 
gravitational potential, still the uni- or bi-directionality will someday be a strong signature 
for an extragalactic WIMP population. Additionally, the essentially monochromatic spec- 
trum of extragalactic WIMPs implies (if they are not too heavy) a nuclear recoil spectrum 
that is flat up to the maximum recoil energy. 

For a critical density of extragalactic WIMPs, the count rate in a detector is roughly the 
same as that of fast galactic WIMPs, but their signature is quite different. An extragalactic 
population of particular interest would be WIMPs associated with the Local Group; these 
may have a 10-30 times higher count rate than galactic WIMPs of the same speed, for an 
energy resolution of 1% in the detector. Detection of extragalactic WIMPs would bring 
valuable new insight into the processes of galaxy formation. 

We also wish to point out that the annual modulation of the signal due to the Earth 
moving around the Sun decreases as one goes to higher velocity WIMPs. The Earth's annual 
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motion of 30 km/s is 5% of the escape velocity and hence 5% of the velocity of extragalactic 
WIMPs. 
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APPENDIX A: TRAJECTORIES OF EXTRAGALACTIC WIMPS. 



In this appendix, we mention some details of the calculation of the trajectories of extra- 
galactic WIMPs entering the galaxy. In particular, if a WIMP at infinite distance from the 
galaxy is approaching with one angle a^, we calculate the angle a at which it arrives at the 
Sun. The deviation of the angle from its initial value can be huge. One might assume that 
radial infall through the location of the Sun to the center of the galaxy is a good approxima- 
tion to the motion; however, this assumption is not always true, depending on the galactic 
potential. In an extreme case, the incoming particle can instead have a purely tangential 
velocity at the location of the Sun, with no radial component (towards the Galactic Center) 
at all. We considered only galactic potentials which are spherically symmetric. Specifically, 
we present our results for (i) a Keplerian (point mass) potential, (ii) a Navarro- Frenk- White 
(NFW) potential, and (iii) a logarithmic potential. For comparison, we also consider (iv) 
the free particle case. 

Gravitational potentials: 

(i) As a simple check, we consider a Keplerian potential with 0Kc P ( r ) = —GM/r, a point 
mass at the center of the galaxy. This is not a realistic potential. 

(ii) An NFW potential corresponds to an NFW density profile p = p s (r/r s ) -1 (l+r/r s )~ 2 , 
where r s is a length scale and p s is four times the density at r = r s . Explicitly, 

0NFw(r) = -4nG P y s ln{1+ . r/rS \ (Al) 

r/r s 

We fix the scale r s by imposing that the escape velocity from the galaxy at the Sun position 
r is f csc , a given quantity. This gives the following equation for u s = r /r s , 

In(l + «.) = 1-M (A2) 



1 + U s ' V, 



2 

esc 



Here v c is the galactic circular velocity at the Sun position. Taking v c = 220 km/s and 
f esc = 600 km/s gives u s = 2, in other words r s = r /u s = r /2.00 or r s = 4 kpc if r = 8.5 
kpc. [We are aware that this value of r s is an approximation to what may arrive in a realistic 
scenario.] 

(iii) A logarithmic potential is interesting in this context because it is a simple function 
with constant circular velocity. However, it increases without limit as r — > oo, and so it 
needs to be cut off at a certain radius R to allow for unbound orbits. Hence we use the 
truncated logarithmic potential 



0Lo g (r) 



,R 



-v 2 r —, for r > R , 

r 



, for r < R , 

(A3) 



fixing the cut-off radius R by again imposing that the escape speed at the Sun position is 
f esc . This gives the condition 

UR= T j i =eM^-vlJ{2vl)]. (A4) 
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With f esc and v c as above, ur = 0.0659. 

Arrival directions: Since we assume spherical symmetry for the galactic potential, angular 
momentum is conserved and the orbit is planar. The particle trajectories that pass by the 
Sun lie on a single plane, which is defined by the particle velocity vector at infinity (the 
stream velocity) and the vector from the galactic center to the Sun ro- [The only 
exception is if the WIMPs come towards the Sun from directly behind the galactic center, 
as discussed below.] 

In the plane of the orbit, we define the angles and a. To characterize the direction 
from which the particle entered the Galaxy, we define the angle to be the angle between 
the following two vectors that originate at the Galactic Center: the particle velocity vector 
at infinity towards the Galactic Center, and a vector pointing from the Sun towards the 
Galactic Center (but translated to the location of the Galactic Center). To characterize the 
direction the particle is moving once it reaches the Sun, we define a as the angle between 
the following two vectors that originate at the Sun: the tangent vector of the particle's 
motion as it passes the Sun, and the vector from the Sun towards the Galactic Center. As 
an example, see Figure 4. In order to obtain as a function of a, we need to solve the 
orbit equations. 

Introduce polar coordinates (r, 6) in the plane of the orbit with the origin at the Galactic 
Center. Let r be the distance of the Sun from the Galactic Center, and v ^ be the particle 
velocity at infinity. Then conservation of energy gives 

\r 2 + \rH 2 + 0(r) = (A5) 

where we have chosen the zero of the potential <j)(r) at infinity. Conservation of angular 
momentum gives 

r 2 9 = v r sin a (A6) 



where v = y^l + t> 2 sc is the particle speed at the Sun location r and t> csc = y— 20(r o ) 
is the escape speed at the Sun location. Substituting eq. (|A6|) into eq. (|A5|), and defining 
u = r /r, we have 

du 1 



d9 sin a 



K(u) - u 2 sin 2 a, ( A7) 



where the upper (positive) sign in eq. (|A7|) applies to particles going in, and the lower 
(negative) sign to particles going out. Here 

K{u) = V - ~ 2 f° M (A8) 

is the ratio of the particle kinetic energies at galactocentric distances r = r /u and r . We 
evaluate K{u) for an NFW, a logarithmic, and a Keplerian gravitational potential, and for 
the free particle case. Then we integrate eq. ([AT]) from # = 0to# = 7r — ctoo with initial 
condition w(0) = to find 

J 7r + 6(a) — ^(a), forcosa>0, , 
a °° ~ 1 7T -9(a), forcos«<0, [ } 
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where 



0(a) = sin a/ , (A10) 

" / ° \JK(u) - -u 2 sin 2 o; 

and 

/•»« rf-u 

#(a) = 2sina/ . (All) 

Jo \JK{u) - -u 2 sin 2 o; 

Here u a is the solution of 

K(u a ) -u 2 a sin 2 a = 0; (A12) 

^(a) is the angle between asymptotes for particle trajectories passing at the Sun position 
with speed v . 

We now present the results of calculating function of a. Figure 3 plots the 

four functions obtained for the NFW potential, the logarithmic potential, the Keplerian 
potential, and the free particle case. The axes have been labeled in degrees. 

The graph in Figure 3 can be inverted to obtain the particle direction a at the Sun posi- 
tion once its direction at infinity is given. This inverse function is in general multi- valued: 
for instance, = 20° in Figure 3 corresponds to two values of a for each gravitational 
potential plotted. The physical interpretation is that for these incoming directions of the 
particle stream, particles reach the Sun position from two separate directions: one beam has 
passed to the "left" of the galactic center, the other to the "right". Notice that for some 
stream directions there is only one possible direction for the particle flow at the Sun. It must 
be remembered that except for = the particle directions at the Sun position lie in the 
plane defined by the directions of the galactic center and of the stream velocity at infinity. 
The case ctoo = is special, in that it corresponds to the Sun being on the downstream 
axis. The values of a read from the graph in Figure 3 then give the half-angle of the cone 
of directions from which the particles approach the Sun. 

Particles approaching the galaxy from the direction of M31 with a speed = 70 km/s 
have v = 604.1 km/s and 

aoo = 58?5, ttKcp = {H3?4, 305?0}, a NFW = 88?3, a Log = 81?4. (A13) 

There is only one particle flow for the NFW and the logarithmic potential. In the NFW 
case the particle velocity at the Sun position is almost tangential. This is because in those 
trajectories that pass by the Sun, the Sun is almost at the point of closest approach to the 
galactic center. At this point, the particle radial velocity vanishes, and the particle velocity 
is purely tangential. 

In Figures 4 and 5, we present examples of two orbits for an NFW potential. Figure 4 is 
a typical case with a relatively small deviation of the angle from its value at infinity, while 
Figure 5 is close to the extreme case in which the orbit has a purely tangential component 
at the location of the Sun. 
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FIG. 1. Expected recoil spectra for Milky Way and Local Group WIMPs. Detector and WIMP 
parameters as indicated. The Local Group spectrum corresponds to an angle ip = 155° between the 
velocity of the Local Group WIMPs and the velocity of galactic rotation at the solar neighborhood 
(both velocity vectors starting at the location of the Earth). The Milky Way spectrum reflects a 
velocity distribution truncated at the escape speed (curve marked v esc = 600). The dotted line, 
marked v esc = oo, does not include this truncation. Thick lines are the sum of the galactic (with 
v esc = 600 km/s) and extragalactic spectra. 
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FIG. 2. Same as Figure 1, but for an angle ip = 90° between the velocity of the Local Group 
WIMPs and the velocity of galactic rotation at the solar neighborhood. 
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FIG. 3. Graph relating the particle direction at infinity Oqo to the particle direction at the 
Sun position a. Clearly for any given incoming velocity from infinity, at most two directions in 
a detector are possible. Solid lines are for the NFW potential, dashed lines for the logarithmic 
potential, dash-dotted lines for the Keplerian potential, and dotted lines for the free particle case. 
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FIG. 4. An example of an NFW trajectory with incoming angle at infinity = 225°. The 
lower dot indicates the position of the Sun, the upper dot denotes the galactic center. The arrival 
direction and the incoming direction at infinity are shown by dashed lines. Also shown are the 
angles and a. 
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FIG. 5. The unique NFW trajectory passing by the Sun for a particle wind coming from an 
angle of 58? 5 from the galactic center (i.e. from the direction of the Andromeda galaxy). Notice 
the almost tangential character of the trajectory near the Sun. 
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